Liquid-liquid phase transition model incorporating evidence for ferroelectric state 
near the lambda-point anomaly in supercooled water. 



P.O. Fedichev 

Quantum Pharmaceuticals Ltd, III. Kosmonavta Volkova 6a-1205, 
125171, Moscow, Russian Federation, \http://q-pharm.~co n "*] 



L.I. Menshikov 

NRC Kurchatov Institute, Kurchatov Square 1, 123182, Moscow, Russian Federation 

We propose a unified model combining the first-order liquid-liquid and the second-order ferroelec- 
tric phase transitions models and explaining various features of the A-point of liquid water within 
a single theoretical framework. It becomes clear within the proposed model that not only does the 
long-range dipole-dipole interaction of water molecules yield a large value of dielectric constant e at 
room temperatures, our analysis shows that the large dipole moment of the water molecules also 
leads to a ferroelectric phase transition at a temperature close to the A— point. Our more refined 
model suggests that the phase transition occurs only in the low density component of the liquid 
and is the origin of the singularity of the dielectric constant recently observed in experiments with 
supercooled liquid water at temperature T ~ 233K. This combined model agrees well with nearly 
every available set of experiments and explains most of the well-known and even recently obtained 
results of MD simulations. 



Supercooled water exhibits a number of well-known 
anomalies near the so-called A— point at temperature 
T\ w 22%K (at normal pressure), where a good number 
of thermodynamic quantities such as heat capacity, com- 
pressibility, thermal expansivity, and dielectric constant 
all exhibit nearly singular behavior [TH3]. The intrinsic 
thermodynamic instability of liquid water at tempera- 
tures well below the freezing point has been a major ob- 
stacle both in experimental studies and theoretical mod- 
eling (see e.g. [1HTT] for a review). The observed features 
are often weak, which suggests a thermodynamic conti- 
nuity of the various water states near the A— point and 
relates the observed "singularity" with essentially a ran- 
dom sharp feature, a function of the parameters of the 
liquid |12l [T5] . A considerably more accepted view is to 
ascribe the features characteristic of a phase-transition to 
a first-order liquid- liquid phase transition (LLPT). This 
model predicts the existence of a second critical point 
of water at the temperature Ten ~ 200K and pressure 
~ lkbar |14l 115] . The sharp temperature dependencies 
observed near T\ are attributed to crossing the Widom 
line [T6], where the density and entropy fluctuations are 
large [17] and which happens at T w T\ at normal pres- 
sure. The view is supported by numerous molecular dy- 
namics (MD) simulations based on realistic water models 
[141 [T8H21)] . simplified analytical models [27-31 , and ex- 
perimental studies |H I32H34] . The reported anomalies 
are not restricted to static features, the dynamic proper- 
ties such as the Einstein relation between diffusion and 
mobility coefficients [3"S"1 l3"6"] and the Arrhenius behavior 
of the liquid's dynamic properties [16, 36 37 break down 
near T\ as well. 

Recent MD studies have demonstrated that heat ca- 
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pacity and thermal conductivity |38| peak around T\ 
as well and that the liquid shows a good deal of or- 
dering in the vicinity of the A— point. Similarly, re- 
cent measurements of the dielectric properties of liquid 
water confined in nanopores and hence prevented from 
freezing well below the natural freezing point [39 mani- 
fest a profound bump in the dielectric constant near the 
A— point [40] 141] . These observations paint a richer pic- 
ture than a mere first-order liquid-liquid phase transi- 
tion and, in fact, bring back an old idea [3] relating the 
weak singularity of the dielectric constant to a ferroelec- 
tric phase transition (FPT). Remarkably, the hypothesis 
was put forward immediately after the discovery of the 
A— point, though the weakness of the observed singular- 
ities prompted the authors [3] to reject the explanation. 
Furthermore, the ferroelectric instability for a model liq- 
uid with parameters similar to water is predicted to occur 
at a very high and essentially unreachable temperature, 
— 1200-ftf A more sophisticated model [4"3"ti4"5] pre- 
dicts a ferroelectric phase transition at a temperature 
independent of the details of the short-range interactions 
between the molecules, 
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where hq is the density and do are the static dipole mo- 
ments of the molecule comprising the liquid. The dielec- 
tric constant is not associated with the molecules' de- 
grees of orientational freedom and comes from electron 
shell polarization, = A~ 5.5 (as discussed in e.g. [464 - 
I48|). Despite bringing the phase transition temperature 
into the right range, the model predicts a behavior of 
the dielectric constant with changes in temperature that 
is far too gradual when compared with empirical results 
[401 |4T] . Consequently, it cannot even qualitatively ex- 
plain all of the features of supercooled water by itself. 
Realistic liquids such as water are far more complicated 
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than a model polar liquid consisting of point dipoles. For 
example, the tetrahedral geometry of H2O molecule and 
its electron shell leads to a polyamorphism phenomenon 
[U [7J US] , namely two or more phases of the same liquid 
existing in a mixture at the same time. 

To unify the observed ferroelectric-like properties of 
water- molecule ordering and the singularity of the dielec- 
tric constant near T\ |40II41) with the previously reported 
signatures of the LLPT within a single theoretical frame- 
work, we combine our simple polar liquid phenomenology 
gDl SHH] with the LLPT hypothesis H [H H] using a 
two-component mixture model of water |50| . We assume 
that the equilibrium state of supercooled water is a mix- 
ture of macroscopically-sized clusters of the two types: 
low density (LDL) and high density liquid (HDL). The 
LDL local lattice is softer than that of HDL and the den- 
sity of HDL exceeds that of LDL by - 20% H [5TH55]. 
Since LDL is "softer", the molecules of the LDL rotate 
more or less freely, whereas in the HDL the rotations 
are more difficult. This explains why HDL has no ferro- 
electric state at any temperature. We assume that the 
ferroelectric ordering and the FPT apparently observed 
in supercooled liquid bulk water occurs in the LDL com- 
ponent only. 

The Gibbs free energy of an LDL cluster at a given 
pressure P is the sum of the contributions from the polar 
liquid, G^ DL , and the lattice, G\ DL : Gldl = G£ DL + 
G\ DL . At high temperatures, T > Tp, the equilibrium 
state of LDL corresponds to the disordered paraelectric 
phase, whereas at lower temperatures, T < Tp. LDL 
undergoes a second-order phase transition and enters the 
long-range-ordered ferroelectric state. Near the phase 
transition, where r = (T — Tp)/Tp <C 1, the free energy 
of LDL takes the form: 

G ldl ^-Dt 2 8(-t) + G l ldl (T,P), (2) 

where D ~ V^n 2 d 2 ~ 150 cal/mol, and Vo — Na/tiq ~ 
22 cm 3 is the molar volume of LDL. Since the LDL den- 
sity differs by no more than 20% from the total liquid 
density, we will not distinguish between the LDL density 
and no- The second-order FPT in LDL manifests itself 
as a singularity of the dielectric constant e (see e.g. |56|): 

e = e 00 (l + /(T)), (3) 

where / (T) = 3T F /(T - Tp) at T > T F and / (T) = 
3T F /2(Tp - T) at T < T F , which is a much weaker de- 
pendence than that observed experimentally [40J. The 
discrepancy is apparently due to the very sharp temper- 
ature dependence of the LDL fraction c near T ss T\. To 
see that, we follow [50 and formulate a two-liquid model 
representing the Gibbs energy of water as the energy of 
a two-liquid mixture of macroscopic clusters using a rep- 
resentation similar to that used in the physics of binary 
alloys (see e.g. [57]): 

G (c) = cG L dl + (1 - c) Ghdl + Uc(1-c) + 

+ RT [clogc +(l-c) log (1-c)]. (4) 
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Figure 1: Model PT diagram of water (see the explanations 
in the text). 



Here Ghdl is the free energy of the HDL component, 
and the parameter U characterizes the "energy of mix- 
ing". The equilibrium free energy and the LDL frac- 
tion c are found through minimization of G (c) over c, 
G 1 (c) = [8G (c, P, T) I dc\ p rp = 0, or 

AG (P, T) +17(1 — 2c) + RT log (j^A = 0, (5) 

where AG(P,T) = Gldl — Ghdl- The temperature of 
the second critical point in the model is Tqr — U / (2R) 

M- 

To analyze these equations, we follow [SU] E3 an d as- 
sume that the thermodynamic quantities in AG (P, T) = 
AP° - TAS° + PAV° are practically temperature- 
independent: 

AP° (P, T) , AS" (P, T) , AV° (P, T) « const. (6) 

The PT diagram suggested by this model is shown 
in Figjl] The AKB line in the Figure is the spin- 
odal line corresponding to G' (c) = 0, G" (c) = 
[d 2 G(c,P,T) /dc 2 ] pT = 0. The KB line is the LDL 
spinodal, where the LDL phase loses its thermodynamic 
stability (the local minimum of the function G (c) corre- 
sponding to the LDL phase disappears on this line). In 
turn, the KA line is the HDL spinodal. The Gibbs po- 
tential G (c) has a single minimum everywhere above the 
AKB line and two minima below this line, at c — c\ and 
c = c 2 . The section KC of the straight line WC is the 
liquid-liquid first-order phase transition line correspond- 
ing to the phase equilibrium conditions: G (ci) = G (02), 
G'(ci) = G'(c 2 ) = 0. According to Eqs.Q and 
these conditions are equivalent to AG (P,T) = 0, which, 
according to approximation (JsJ , means that T is a linear 
function of P on the KC line, which means that the KC 
line is itself a straight line in Figure [l] The low and high 
density amorphous ice regions are denoted by LDA and 
HDA. 
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Figure 2: Calculated water heat capacity (formulas Q, (TsJ) , 
(10 1; solid curve) versus MD calculations |38| (points) and 
empirical data 1 (crosses). 



The KD line in Figure [T] corresponds to the limit 
AS — > at T — > as required by the third law of ther- 
modynamics and which is never the case under the as- 
sumptions pj. WK is the Widom line, defined by the 
conditions G7 (c) = and G'" (c) = 0, or, equivalently. 
c = 1/2. At small pressures the Widom temperature is 
T w — AE°/AS°. The heat capacity, Cp, consists of the 
two parts 
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where Cp is the lattice contribution 



calculated e.g. 



which can be 
using Debye approximation 
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= 18cal ■ 

mol K (as suggested in [T]) and 9 « 150K is the 
Debye frequency of water. Within the described model 
the polar part of Cp is 
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At small pressures along the Widom line c sa 1/2, G" (c) 
is small, the fluctuations are strong, and the temperature 
dependence of the heat capacity contribution 
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takes a standard Lorenz form, where A = Tw — Tqr, S — 
2RA V / A/(AE°AS°). The quantity peaks at T = T w 
in agreement with experiments [51 [5] , earlier explanations 
|16| . and recent calculations [38 . Therefore we can use 
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Figure 3: Theoretical temperature dependence of the LDL 
fraction (solid lines, see the explanations in the text) versus 
empirical data [33] (dots). 



Eq. ^ to analyze the heat capacity calculated e.g. using 
MD from [35] and extract the model parameters: 

AE° = -860cal ■ moV 1 , AS = S.bcal ■ mol^R- 1 , 

and U = 880cal ■ mol . These parameters correspond 
to T w = 2A5K, T CR = 220K, and P CR = lkbar. 
To calculate Pcr we used AV° w 3.8cm 3 mol~ 1 from 
|50| . where very similar parameters were obtained: U ~ 
900cal ■ mol' 1 , AE° = -250cal ■ mol' 1 , AS = -leal ■ 
mol^K' 1 , and T CR = 225K. These "original" pa- 
rameters yield smaller critical pressure value: Pcr ~ 
0.33kbar. The heat capacity Cp calculated with the help 
of Eqs. 0, ^ and ^ is plotted on Fig{2] against the 
experimental values from [l j and the recent MD calcula- 
tions [35]. Note that although the MD calculation does 
not provide a full match with the experimental curve 
in absolute terms, both data sets consistently describe 
the same feature and hence apparently the same physics. 
This in combination with the water molecules ordering 
predicted by the model (the entropy of LDL is less than 
that of HDL, A5° < 0, in accordance with [B]) and con- 
firmed by the MD calculations is the indication of ferro- 
electric transition implicitly present already in the simu- 
lation 138 . 



Once the model parameters ( 10 1 are established we 



can verify the consistency of the model by observing the 
temperature dependence of the LDL fraction c (T) given 
by Eq. ([5| (see solid line 1 in Figure |3| and compare it 
with the empirical data [33] (the dots in the same Fig- 
ure). There is a qualitative agreement at least at suffi- 
ciently large temperatures above 200K. We note, that 
AG(P,T = 0) < as T -> and therefore, at suffi- 
ciently low temperatures all of the liquid should turn 
into LDL, i.e. c — >• 1 as T — > 0. This leads us to be- 
lieve that the equilibrium composition of water was not 
actually achieved in the experiments of [3*3"] , even at very 
low temperatures. This may well be due to the increase 
of equilibration time at extremely low temperatures, as 
discussed in e.g. [551 [5l3] . With this in mind, we can at- 
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Figure 4: Dielectric const ant fo r bu lk water versus tempera- 
ture. Results of formulas ||12[), ( |10[ ) are given by solid curve. 
Experimental results 40, 41 are presented by points, crosses 
correspond to the measurements [3J. 



tempt to match the LDL fraction measurements only at 
higher temperatures and obtain better agreement with 
the experiment using slightly different model parameters 

AE° = ~920cal ■ mol' 1 , AS = -Acal ■ mol^K' 1 , 

(11) 

and U — 900cal ■ mol 1 from [50] (see solid line 2 in 
Figure [3]). Although curves 1 and 2 are clearly different, 
both sets, ( 10 ) and (111, are very close numerically, which 



indicates further difficulties in refining the parameters 
due to a very sharp dependence of the LDL fraction c next 
to the Widom line and clear experimental difficulties. 

An analysis of the heat capacity and the LDL frac- 
tion measurement let us verify the model and deter- 
mine the model parameters. Consequently, we may use 
the model to predict the temperature dependence of 
the dielectric constant and compare it with empirical 
data. Accordingly, the dielectric constant of the liq- 
uid is the sum of the LDL, cldl = £oo(l + f(T)) with 
f(T) from Eq. and the HDL contributions. The 

HDL fraction is 1 — c, and its dielectric constant can 
be described by the Debye-Onsager model |60| : chdl = 
eoo + 2im dl ( eoo + 2) /(9T). Therefore the full expres- 
sion e = ctLDL + (1 — c)£hdl becomes 



= £oo [l + c(T)f(r)} + [l-c(T)} 
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The predicted temperature dependence calculated using 
the parameters (JToJ) and the value eoo =4.7 correspond- 
ing to Tp = 233K is compared with recent experimental 
measurements [3] HOI IH] in Figure |4j At high tempera- 
tures, T > Tp, Ty/ the model dependence is not far from 
the experimentally observed values. The temperature 
dependence is now sharp, much stronger than that pre- 
dicted by a simpler, single-component model [40 . This 



means that including more liquid states in the analy- 
sis clearly improves the agreement with experimental re- 
sults. Below the transition point to the left of the peak, 
at T < Tp, experiments yield distinctly different results. 
There may be many reasons for that. For example it is 
quite possible there is more than one LDL state of water 
below Tp. The temperature Tp itself may depend on the 
LDL cluster size, and therefore we may face a "continu- 
ous set" of ferroelectric transitions in a multitude of LDL 
forms as the temperature decreases. 

When compared with available empirical data and nu- 
merical calculations, the model calculations of the heat 
capacity, LDL fraction, and the dielectric constant sup- 
port a broader view, implying that supercooled water is 
indeed a mixture of at least two different components, 
namely LDL and HDL. There is a growing body of ev- 
idence for ferroelectric transition in LDL clusters. In- 
dependent of the theoretical arguments given above, the 
FPT in LDL can be supported by the structural similar- 
ity between the local structure of LDL and the crystalline 
lattice of stable ice Ih [oTl 152"] . At normal pressure the 
paraelectric ice Ih enters the ferroelectric state, ice XI, at 
temperature T — 72K [63-66] . The ferroelectric state, 
ice XI, is obtained from hexagonal ice Ih after reconstruc- 
tion of the crystalline lattice, namely, by the shifting of 
water molecules in each elementary cell of the crystal. It 
is interesting to note that theory |67| predicts antiferro- 
electric ordering of molecular dipole moments for Ih ice 
with undeformed crystalline lattice, but the ferroelectric 
state with a deformed lattice has a lower free energy. It is 
not surprising that, due to a less dense local crystalline 
lattice of LDL, the transition temperature observed in 
|40l |4"T] essentially exceeds the corresponding value of 
FPT "ice Ih^=ice XI". It is also worth mentioning here 
that there have been important studies concerning the 
ferroelectric states of metastable ice Ic, with a cubic lat- 
tice first predicted in |68| and detected experimentally 
in [69} [70]. The reported explanation of the "ice Ih^=ice 
XI" transition gives us cause to believe that the ferro- 
electric phase transition in LDL is also accompanied by 
the reconstruction of the local crystalline lattice, i.e. the 
dipolemoment orientation of molecules is strongly linked 
to the lattice's degrees of freedom. On the other hand, 
FPT occurs only in LDL. Therefore, it is reasonable to 
assume that the KC line in Figure [T] is in fact the bound- 
ary between the ferroelectric and paraelectric regions of 
the supercooled water. Here |AGp| ~ |AG/,| and there- 
fore, the molecular orientation and the lattice's degrees 
of freedom are equally important "driving forces" of both 
the ferroelectric and the liquid-liquid phase transitions. 
Above the FK line the situations is different; there typi- 
cally |t| < 1 andEq. ^ yields |AG P | <C |AG L |. There- 
fore, the lattice degrees of freedom dominate. Hence, we 
believe that supercooled water is ferroelectric below the 
FKP line. 
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